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Summary 

-ru oneri c tf 1 8 team focused its efforts on using remotely sensed data to characterize the 

r& B ,Cw^S^ca. Cassificahon of ,he SSA. The dara are 

Pr0V N® toSSTheLu e sSeso„ rhe BOREAS CD-ROMs have been compressed using the 
Gzip program. See Section 8.2 for details. 


Table of Contents 

1) Data Set Overview 

2) Investigator(s) 

3) Theory of Measurements 

4) Equipment 

5) Data Acquisition Methods 

6) Observations 

7) Data Description 

8) Data Organization 

9) Data Manipulations 

10) Errors 

11) Notes 

1 2) Application of the Data Set 

13) Future Modifications and Plans 

14) Software 

15) Data Access 

16) Output Products and Availability 

17) References 

18) Glossary of Terms 

19) List of Acronyms 

20) Document Information 


1. Data Set Overview 


1.1 Data Set Identification 

BOREAS TE-18 Biomass Density Image of the 88 A 

02-Sep-1994 based on a technique developed by Dr. Forrest Han. 
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1.3 Objective/Purpose 

1.4 Summary of Parameters 

Biomass Density 

1.5 Discussion 

mrhH tanCeS |'? f thrCe Cnd n l embers that make U P the land surface. The three main surface features 
c^odv mino?^’ '“A 11 background ’. and shadow - The points of the trajectory range from 0% 
n r ° A geometric optical canopy model was used to determine the areal 

proportions of each of these elements. The nearest trajectory point to each pixel of the lSL TM 
■mage was de.erm.ned to derive ,he biomass densiiy LedY ihe a« 


1.6 Related Data Sets 

BOREAS TE- 1 8 Landsat TM 
BOREAS TE- 1 8 Landsat TM 
BOREAS TE- 18 Landsat TM 


Maximum Likelihood Classification Image of the NSA 
Maximum Likelihood Classification Image of the SSA 
Physical Classification Image of the SSA 


2. Investigator(s) 

2.1 Investigator(s) Name and Title 
Dr. Forrest Hall 

Biospheric Sciences Branch 

National Aeronautics and Space Administration (NASA) 

Goddard Space Flight Center (GSFC) 

2.2 Title of Investigation 

TE-18 Regional Scale Carbon Flux from Modeling and Remote Sensing 

2.3 Contact Information 

Contact 1: 

Dr. Forrest G. Hall 
NASA GSFC 
Code 923 

Greenbelt, MD 2077 1 
(301)286-2974 
(301)286-0239 (fax) 

Forrest.G.HaIl@gsfc.nasa.gov 
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Contact 2: 

David Knapp 
Raytheon ITSS 
NASA GSFC 
Code 923 

Greenbelt, MD 2077 1 
(301)286-1424 
(301)286-0239 (fax) 

David.Knapp@gsfc.nasa.gov 


3. Theory of Measurements 

This biomass density image was derived from a Landsat-5 TM scene collected on 02-Sep-1994. A 
classification image was produced using end member reflectance trajectories. Please refer to the 
document for the BOREAS Terrestrial Ecology (TE)-18 Landsat TM physical classification image ot 

the BOREAS SSA. J . . ...... ... 

The general theory behind the method for deriving biomass density is very similar to the way that 
the classification image was created. The reflectance of a pixel is represented in this method as a linear 
combination of the reflectance of sunlit canopy, sunlit background, and shadow. A model or the 
canopy geometries was used for the various classes. This model computes the proportions of canopy, 
background, and shadow that occupy the area in a pixel based on solar and canopy geometry Using 
ground and satellite reflectance measurements and the proportions of the end members, model 
reflectances were produced for each class. These model reflectances form a trajectory in 
red/near-infrared space within each class as the proportions of canopy cover go from 0% to lUU/o. 

Ground measurements of biomass density were made for several sites in the SSA. These 
measurements were used to define a relationship between the biomass density and the amount ot 
canopy cover at each site. A scaling factor was determined based on this relationship. The amount ot 
canopy cover at each pixel was determined by finding the model reflectance that has the shortest 
distance in red/near-infrared space. The scaling factor is then applied to the amount of canopy cover tor 
each pixel, resulting in an estimate of biomass density for that pixel. 


4. Equipment 


4.1 Instrument Description 

The Landsat-5 TM sensor system records radiation in seven bands ot the electromagnetic 
spectrum. It has a telescope that directs the incoming radiant flux obtained along a scan line through a 

scan line collector to the visible and near-infrared focal plane, or to the mid-infrared and 

thermal-infrared cooled focal plane. The detectors for the visible and near-infrared bands (1-4) are tour 
staggered linear arrays, each containing 1 6 silicon detectors. The two mid-infrared detectors are 1 6 
indtum-antimonide cells in a staggered linear array, and the thermal-infrared detector is a four-element 
array of mercury-cadmium-telluride cells. 

4.1.1 Collection Environment 

The data used to produce this classification were collected by the Landsat-5 TM on 02-Sep-1994. 
Landsat-5 orbits Earth at an altitude of approximately 705 kilometers. 

4.1.2 Source/Platform 

Landsat-5 satellite 
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4.1.3 Source/Platform Mission Objectives 

The mission of the Landsat-5 satellite is to measure reflected radiation from Earth's surface at a 
spatial resolution of 30 meters and to measure the temperature of Earth's surface at a spatial resolution 
of 120 meters. 

4.1.4 Key Variables 

Reflected radiation Emitted radiation Temperature 

4.1.5 Principles of Operation 

The TM is a scanning optical sensor operating in the visible and infrared wavelengths. It contains a 
scan mirror assembly that directly projects the reflected Earth radiation onto detectors arrayed in two 
focal planes. The TM achieves better imagery resolution, sharper color separation, and greater inflight 
geometric and radiometric accuracy for seven spectral bands simultaneously than the previous 
Multispectral Scanner (MSS). Data collected by the sensor are transmitted to Earth-receiving stations 
for processing. 

4.1.6 Sensor/Instrument Measurement Geometry 

The TM depends on the forward motion of the spacecraft for the along-track scan and uses a 
moving mirror assembly to scan in the cross-track direction (perpendicular to the spacecraft). The 
Instantaneous Field of View (IFOV) for each detector from bands 1 -5 and band 7 is equivalent to a 
30-meter square when projected to the ground; band 6 (the thermal-infrared band) has an IFOV 
equivalent to a 1 20-meter square. 

4.1.7 Manufacturer of Sensor/Instrument 

NASA GSFC 
Greenbelt, MD 2077 1 

Hughes Aircraft Corporation 
Santa Barbara, CA 

4.2 Calibration 

The internal calibrator, a flex-pivot-mounted shutter assembly, is synchronized with the scan 
mirror, oscillating at the same 7-Hz frequency. During the turnaround period of the scan mirror, the 
shutter introduces the calibration source energy and a black direct-current restoration surface into the 
100 detector fields of view. 

The calibration signals for bands 1-5 and band 7 are derived from three regulated tungsten-filament 
lamps. The calibration source for band 6 is a blackbody with three temperature selections, commanded 
from the ground. The method for transmitting radiation to the moving calibration shutter allows the 
tungsten lamps to provide radiation independently and to contribute proportionately to the illumination 
of all detectors. 
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4.2.1 Specifications 

The following spectral bands are collected by 


the TM sensor: 


Channel Wavelength (pm) 


Primary Use 


2 

3 

4 

5 

6 
7 


0.45 - 0.52 

0.52 - 0.60 
0.63 - 0.69 

0.76 - 0.90 
1.55 - 1.72 

10.4 - 12.5 

2.08 - 2.35 


Coastal water mapping, soil vegetation 
differentiation, deciduous /coniferous 
differentiation. 

Green reflectance by healthy vegetation. 
Chlorophyll absorption for plant species 
differentiation. 

Biomass surveys, water body delineation. 
Vegetation moisture measurement, snow 
cloud differentiation. 

Plant heat stress measurement, other 
thermal mapping. 

Hydrothermal mapping . 


Band 


Radiometric 
Sensitivity [NE (dP ) ] * 


1 

2 

3 

4 

5 

6 
7 


0 . 8 % 

0 . 5 % 

0 . 5 % 

0.5% 

1 . 0 % 

0.5 K [NE(dT) ] 
2 .4% 


Ground IFOV 

Avg. altitude 
Data rate 

Quantization levels 
Orbit angle 
Orbital nodal period 
Scan width 
Scan angle 
Image overlap 


30 m (bands 1-5, 7) 
120 m (band 6) 

699.6 km 
85 Mbps 
256 

8.15 degrees 
98.88 minutes 
185 km 

14.9 degrees 
7.6% 


* N R The radiometric sensitivities are the noise-equivalent reflectance differences for the reflective 

[NE(dP)l and ,=%« differences for .he , henna.-, nfrared 

bands [NE(dT)]. 


The TM channels were designed for a noise-equivalent differential represented by the radiometric 
sensitivity shown in Section 4.2. 1 . 
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4.2.3 Other Calibration Information 

Relative withm-band radiometric calibration, to reduce "striping," is provided by a scene-based 
procedure called histogram equalization. The absolute accuracy and relative precision of this calibration 
thS 6 aS | Um . e f that any change ,n the optics of the primary telescope or the "effective radiance" from 
It rt " ema Cdllbrat or lamps is insignificant in comparison to the changes in detector sensitivity and 
ectronic gam and bias with time and that the scene-dependent sampling is sufficiently precise for the 
required withm-scan destnpmg from histogram equalization. Each TM reflective band^nd the internal 
ca hbrator lamps were calibrated prior to launch using lamps in integrating spheres that were in turn 
calibrated against lamps traceable to calibrated National Bureau of Standards lamps Sometimes the 
absolute radiometric calibration constants in the "short-term" and "long-term parameters" files used for 

rh^nae^ 0 ^ 88111 ^ bavebeen ^odified after launch because of inconsistency within or between bands 
anges in the inherent dynamic range ot the sensors, or a desire to make quantized and calibrated 
values from one sensor match those from another. 4 canorated 


5. Data Acquisition Methods 


Rem«e Sen'll SS? (T "* “T and received from ^ Canadian Cenire for 

tFOS ATI A 8 CC ^ S f ) ’ h ”P ar c chased them from the Earth Observation Satellite Company 
(EOS A T). As received from CCRS, the image had been processed from raw telemetry to a ' 
systematically corrected product within the CCRS MOSAICS system y 


6. Observations 


6.1 Data Notes 

T . m I hi , S r a fa ry a Wa D T lleCted c° n 02_Se P* 1 994. This scene is Path 37, Row 22-23 (shifted) in the 
Landsa Worldwide Reference System (WRS). The solar elevation angle at the time of image 

7,16 S °' ar aZimU ' h WaS 146 ' degrees ' 71,6 radi °™" fc of 

Th f ™ ir [ ia S e / rom which this classification was produced was atmospherically corrected usine 
aerosol optical thickness data measured by sunphotometers in the study area. These^optical thicknes § s 
da a were used in the Second Simulation of the Satellite Signal in the Solar Spectrum ( 6 S) program to 
The^f .m e r 8phenCa alb ^ do ’ P ath fiance, gaseous transmission, and scattering transmission 
Markham eTS o 992") " ° detenmme SUlface reflectance based on equations 4 a and 4 b of 


6.2 Field Notes 

Not applicable. 


7. Data Description 


7.1 Spatial Characteristics 


7.1.1 Spatial Coverage 

• | T |? e Image area that was classified covers an area that is approximately 144 km by 1 14 km and 

BOR a S PHH n H , nnce Albe K rt ’ A Saskatchewan - The comers of the data senate as follows. 
Sm.on 7 E L4 S ° d COOrdinates are in the Albers Equal-Area Conic (AEAC) projection described in 
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BOREAS Grid NAD 8 3 


Corner 

X 

Y 

Long . 

Lat . 

Northwest 

Northeast 

Southwest 

Southeast 

297 . 810 

441 . 810 

297 . 810 

441 . 810 

392 .490 

392.490 

278.490 
278 .490 

106. 401W 
104 . 190W 
106.515W 
104 . 357W 

54 . 438N 
54 . 333N 
53.417N 
53.314N 


7.1.2 Spatial Coverage Map 

Not available. 

7.1.3 Spatial Resolution 

Each pixel represents a 30-meter by 30-meter area on the ground. 


are^mapped is projected in the BOREAS Grid projection, which is based on the ellipsoidal 
version of the AEAC projection. The projection has the following parameters: 


Datum: North American Datum of 1983 (NAD83) 

Ellipsoid: Geodetic Reference System of 1980 

System of 1984 (WGS84 ) 

Origin : 111.000°W 51.000°N 

Standard Parallels: 52° 30* 00"N 

58° 30 ' 00"N 

Units of Measure: kilometers 


(GRS80 ) or Worldwide Geodetic 


7.1.5 Grid Description n , „ 

The data are referenced to the projection described in Section IAA. 

7.2 Temporal Characteristics 

This Imagery w as° collided on 02-Sep-1994. This scene is Path 37, Row 22-23 i (shifted) ir i the 
Landsat WRS. The solar elevation angle at the time of image acquisition was 40.1 degrees. The so 
azimuth angle was 146 degrees. The radiometric quality of this imagery was acceptable. 

7.2.2 Temporal Coverage Map 

Not applicable. 

7.2.3 Temporal Resolution . . . AO c 1QQ4 

This data set represents the biomass density as it existed on 02-Sep-19V4. 

7.3 Data Characteristics 


7.3.1 Parameter/Variable 

Biomass density. 

7.3.2 Variable Description/Definition 

Biomass density: The amount of biomass per unit area. 

7.3.3 Unit of Measurement 

Tenths of kilograms of biomass per square meter. 
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7.3.4 Data Source 

Landsat-5 TM scene on 02-Sep-1994 from the CCRS. 

7.3.5 Data Range 

Data are stored as 8-bit integers ranging from 0 to 255. Values of 254 and 255 represent pixels for 
which biomass density was not computed (i.e., not conifer). Values of 254 represent water pixels 
Valid biomass density values in this image range from 0 to 253 (0.0 to 25.3 kg of overstorv biomass 
per square meter). y 

7.4 Sample Data Record 

Not applicable for image data. 


8. Data Organization 

8.1 Data Granularity 

The smallest amount of data that can be ordered is the entire data set. 

8.2 Data Format 

8.2.1 Uncompressed Data Files 

The SSA classification product contains two files as follows: 

File 1: (80-byte American Standard Code for Information Interchange (ASCII) text records) 

Text file listing the files on tape 

File 2: (3,800 records of 4,800 bytes each) (1 byte per pixel) 

Image with values from 0 to 255 

8.2.2 Compressed CD-ROM Files 

On the BOREAS CD-ROMs, file 1 listed above is stored as ASCII text; however, file 2 has been 
compressed with the Gzip compression program (file name *.gz). These data have been compressed 
usmg gzip version 1.2.4 and the high compression (-9) option (Copyright (C) 1992-1993 Jean-loup 
Gailly). Gzip (GNU zip) uses the Lempel-Ziv algorithm (Welch, 1994) used in the zip and PKZIP 
programs The compressed files may be uncompressed using gzip (-d option) or gunzip. Gzip is 
available from many Web sites (for example, ftp site prep.ai.mit.edu/pub/gnu/gzip-*. *) for a variety of 
operating systems in both executable and source code form. Versions of the decompression software 
tor various systems are included on the CD-ROMs. 


9. Data Manipulations 

9.1 Formulae 

Not applicable. 

9.1.1 Derivation Techniques and Algorithms 

The techniques used to create this biomass image are described in Sections 1 .5, 3, and 6. 1 

9.2 Data Processing Sequence 
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’V P ^e TM?mage™ P wa S converted to surface reflectance before the Womaas densjty irn^e was 

created* Atmospheric correction coefficients were computed with 6S ustng opt, cal depths from 
a sunphotometer (Markham et al., 1992). 

Tr^^^i^'wer^comfm^e^ba^d on end member reflectances, solar geometry, tree height to 

. Grom^meksIIremenK ™ biomass density were used to find the relationship between biomass 

. The^raiectories^nd^hetmage classification were used as input to the image biomass program. 
I Se bkJmassTmage was mapped into the ABAC projection using nearest ne.ghbor resamphng. 

I staff copied the ASCI and compressed the bina^ fiie 

for release on CD-ROM. 


9.2.2 Processing Changes 
None. 


9.3 Calculations 

9.3.1 Special Corrections/Adjustments 

None. 

9.3.2 Calculated Variables 

None. 

9.4 Graphs and Plots 

None. 


10. Errors 


of various features that tall within a 30 -meter pixel. 

10.2 Qualit} Assessment 

I0 'The biiam-' w^snM checked "(various locations, and the image class was compared to the forest 

= P s^ 

SresemsTteo nm^prerenleTby sites, 

to be someWusion ta tween .claves 

Sy CS S a ^ 

made m minimise this error. it can still occur in pixels .ha, do no, match the model output. 
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10*2.3 Measurement Error for Parameters 

_ Co ™P ans on ^ ith ground measurements of biomass density showed that the root mean sou a red 
error of biomass density for this image is approximately 3.75 kg/m2 where the tree species is very 
homogenous However, errors can vary significantly because of variations in pixel reflectance This 
techmque is thought to be reasonably accurate, especially in areas where blackspmce is relatively pure 

“ Signing W °' her speaes ' areas where different ,ree S P^ ™ogeto eS c^n 

10.2.4 Additional Quality Assessments 

None. 

10.2.5 Data Verification by Data Center 

reasonable 3861 ^ W “ ViSUa " y CheCked &t Vari ° US loCations l ° See that the biomass density values were 


11. Notes 

11.1 Limitations of the Data 

This data set is based on an image that was collected on 02-Sep- 1 994 and represents the biomass 

see Sec,ion l0 - 23 “ 

11.2 Known Problems With the Data 

water'ehss classification show U P in the disturbed class, and cloud shadows show up in the 

ater class. The scene is mostly clear, so this problem has a very limited impact. P 

11.3 Usage Guidance 

ho.H B f 0re UnCOmpre ' S j‘" g the - Gzip flles on CD-ROM, be sure that you have enough disk space to 
CD-ROmZ ySm Then ” ,hC aPPr ° Pria ' e dcc ° m P"*"°" P™*™ Provided on .he 

11.4 Other Relevant Information 

None. 


12. Application of the Data Set 

BOREAS SSA U * d ^ m ° deling PUrp0SeS ‘° ge ' “ r ° Ugh idea of ,he variation in biomass d =nsity 

13. Future Modifications and Plans 

None. 
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14. Software 


14.1 Software Description under EASI/PACE image processing software from PCI, 

UmpS Ziv algorithm (Welch, 1994) used in the zip and PKZIP commands. 

EASI/P A Cdi i^a proprietary software package developed by PCI, Inc. Comae, PCI for details. 


PCI, Inc. 

50 West Wilmot St. 

Richmond Hill Ontario, Canada L4B 1 M5 

(905)764-0614 

(905) 764-9604 (fax) 




CD-ROMs. 


15. Data Access 


The SSA biomass density image data set is 
Information System (EOSDIS) Oak Ridge 
Center (DA AC). 


ivailable from the Earth Observing System Data and 
lational Laboratory (ORNL) Distributed Active Archive 


15.1 Contact Information 

For BOREAS data and documentation please contact. 


ORNL DAAC User Services 
Oak Ridge National Laboratory 
P.O. Box 2008 MS-6407 
Oak Ridge, TN 37831-6407 
Phone: (423) 241-3952 
Fax: (423) 574-4665 

E-mail: omldaac@oml.gov or oml@eos.nasa.gov 


,5 -L D rih o£"!ng System Damand Information System (EOSDIS) Oak Rtdge National Laboratory 
(ORNL)* Distributed Active Archive Center (DAAC) for Btogeochemtcal Dynanucs 

http://www-eosdis.oml.gov/. 

* U^rs may” abtain°data dire'ctly'through the ORNL. ^ or by 

mSKephone, fax, letter, or personal visit using the contact 

information in Section 15.1. 
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15.4 Data Center Status/Plans 

ha^opy°?am L p^u A cR 'Ue goREASCD ROM® S55“ fle f ld "“T T™' ™age, GIS, and 
CD-ROM are available tarn DAA? d “ ° T l,S,ed in inven,ories the 


16 . Output Products and Availability 

16.1 Tape Products 

or 625ol*es fS tachW^ 1 *^ °" S ' mm ' Disi,a ' Archive Tape <DAT) - or 9 - ,rack “P“ « 1600 

16.2 Film Products 

None. 


16.3 Other Products 

These data are available on the BOREAS CD-ROM series. 


i /. Keterences 

H 7 ail p^ tf ® r ™^ ensor/I " s ‘ r “ment/Data Processing Documentation 

PACE Image Analysis Kemal Version 5.2. 1993. PCI, Inc. Richmond Hill, Ontario. 

Richards, J.A. 1986. Remote Sensing Digital Image Analysis: An Introduction. Springer Verlag. 
No'^VM-I? 4 ’ A TCChniqUe f ° r High Per f orrnan ce Data Compression. IEEE Computer, Vol. 17, 

1^*2 Journal Articles and Study Reports 

I WMA’naSA KR5'(&m 0SPh ™ S,Udy: Experimen ‘ Plan ' Versi0n 

ss B^^ R B e^' ( r x 'ffi 6 r osphere s,udy: Expenmc "' pian ' v -» 
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Sellers, P„ F. Hall, and K.F. Huenrmnch M ^Boreal EcosyslenvABrosphere Study: ,994 
Operations. NASA BOREAS Report (OPS DOC 94). 

Sellers, P„ F. Hall, and K.F. HuenrmnchM 997. Boreal Ecosystem-Atmosphere Study: 1996 
Operations. NASA BOREAS Report (OPS DOC 96). 

Sellers, P., F. Hall, H. Margolis, B. Kelly, D Jaldocchu a den Hartog, J. Cmiar, M.O. Ryan, B. 
“ ^ 1994 fie,d year - 
Bulletin of the American Meteorological Society. 76(9). 1 549 1 S • 

Sellers, P.J., F.O. Hall, R.D Kelly. A Black D. Baldocclu ^ 

Crill, D.P. Lettenmaier, H. Margolis, J Cihlar ). Ne • J p g Guertin 1997 . BOREAS 

in° 1 ^ 7 :^x^ri^n*t (^e^i ld w^cientific Results arid Future Directions. Journal of Geophysical 
Research 102 (D24): 28,731-28,770. 

17.3 Archive/DBMS Usage Documentation 

None. 

18. Glossary of Terms 

None. 


6S 

AEAC 

ASCII 

BOREAS 

BORIS 

BP I 

OCRS 

CD-ROM 

DAAC 

DAT 

DEM 

EOS 

EOSAT 

EOSDIS 

CIS 

GMT 

GRS80 

GSFC 

IFOV 

MSA 

MSS 

NAD2 7 

NAD 8 3 

NASA 

NSA 

ORNB 


19. List of Acronyms 

... <-u„ qatellite Signal in the Solar Spectrum 

■ second Simulation of the Satellite oiy 

• Albers Equal-Area Conic 

. American standard Code for Information Interchange 

- BOReal Ecosystem-Atmosphere Study 

- BOREAS Information System 

- Bytes Per Inch 

- Canadian Centre for Remote Sensing 

- Compact Disk-Read-Only Memory 

- Distributed Active Archive Center 

- Digital Archive Tape 

- Digital Elevation Model 

- Earth Observing System 

- Earth Observing Satellite Company 

- EOS Data and Information System 

- Geographic Information System 

- Greenwich Mean Time 

- Geodetic Reference System of 1980 

- Goddard Space Flight Center 

- instantaneous Field of View 

- Modeling Sub-Area 

- Multispectral Scanner 

- North American Datum of 1927 

- North American Datum of 1983 

- National Aeronautics and Space Administration 

- Northern Study Area 

- Oak Ridge National Laboratory 
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PANP 

RSS 

SSA 

TE 

TM 

URL 

UTM 

WGS84 

WRS 

WWW 


- Prince Albert National Park 
Remote Sensing Science 
Southern Study Area 

- Terrestrial Ecology 

- Thematic Mapper 

- Uniform Resource Locator 
Universal Transverse Mercator 

- World Geodetic System of 1984 
Worldwide Reference System 

- World Wide Web 
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